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ABSTRACT 

The known host galaxies of short-hard gamma-ray bursts (GRBs) to date are characterized by low 
to moderate star-formation rates and a broad range of stellar masses. In this paper, we positionally 
associate the recent unambiguously short-hard Swift GRB 100206A with a disk galaxy at redshift 
z = 0.4068 that is rapidly forming stars at a rate of ~ 30 Mq yr~^, almost an order of magnitude higher 
than any previously identified short GRB host. Using photometry from Gemini, Keck, PAIRITEL, 
and WISE, we show that the galaxy is very red {g — K = 4.3 AB mag), heavily obscured {Ay ~ 2 mag), 
and has the highest metallicity of any GRB host to date (12 -I- log[O/H]KD02 = 9.2): it is a classical 
luminous infrared galaxy (LIRG), with iiR w 4 x lO-^^L©. While these properties could be interpreted 
to support an association of this GRB with very recent star formation, modeling of the broadband 
spectral energy distribution also indicates that a substantial stellar mass of mostly older stars is 
present. The current specific star-formation rate is modest (specific SFR « 0.5 Gyr~^), the current 
star-formation rate is not substantially elevated above its long-term average, and the host morphology 
shows no sign of recent merger activity. Our observations are therefore equally consistent with an older 
progenitor, similar to what is inferred for other short-hard GRBs. Given the precedent established 
by previous short GRB hosts and the significant fraction of the Universe's stellar mass in LIRG-like 
systems at z > 0.3, an older progenitor represents the most likely origin of this event. 
Subject headings: gamma-ray burst: individual: 100206A 



1. INTRODUCTION 

Classical gamma-ray bursts (GRBs) segregate into 
two phenomenological classes on the basis of their 
prompt-emission properties: long-duration, soft- 
spectru m GRBs and short-d uration, hard-spectrum 
GRBs (|Kouveliotou et al.l[T993[ ). The two distributions 
overlap in both duration and in hardness, but a duration 
of Tgo = 2 s is commonly taken as the dividing line 
between "long" and "short" bursts^. Long bursts 
constitute the large majority of events detected by all 
major GRB sateUites (about 75% of BATSE events and 
over 95% of Swift events); because these events hav e 
brighter and longer- lived afterglows (|Kann et al.ll2011[ ). 
they have historically been easier to study. Observations 
stretching back almost 15 years associate long-duration 
bursts exclusively with actively star-forming host galax- 
ies, and specifically with regions bright at ultraviolet 
(UV) wavelengths (Fruchter et all l2006f ) close to the 
center (Hi oom et al.l l2002f ) of these galaxies, strongly 
pointing toward a massive stellar origin. In a number 
of cases, evidence for an accompanying supernova has 
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However, the question of classification has recently become 
much more com plex f ollowing the discovery of events with ex- 
tended emission l INorris fc B onncll 2006) and deep limits on supcr- 



provided direct ev idence confirming this conclusion (see 
iWooslev fc Bloo^ [2b06 for a review . 

The short GRB sample is much smaller, but its 
properties are unambiguously different from those 
of long bursts. At least two short events (GRBs 
050509 B and 050724: IBloom et al.|[200l IGehrels et al.l 
2005t iBarthelmy et all 120051: Berger et al.l 120051 : 
Prochaska et al. I2006D have been associated with 



IGal-Yam eralll200g ). bee 

the subject of classification 



novae from t > 2 s b ursts i Fynbo ct al. 2006; Gehrcls ct al. 2006; 



Bloom et al., (,2008i ) for a discourse on 



massive, evolved galaxies having essentially negligible 
current star formatio n, pointing strongly tow ard a 
long-lived progenitor (IBloom fc Prochaskal l200"6l ). In 
addition, the offset distribution for short GRBs relative 
to their host galaxies extends to much greater distances 
than for long GRBs — partially as a simple result of the 
fact tha t short-GRB host s tend to be physically more ex- 
tended ()Fong et al.ll2010D , but short GRBs also generally 
show minimal correlation with the host blue /UV light 
or re s olved regions of active star format i on (|Fox et al.l 
I2005t iFong et al.l l2010t iRowlinson eFall l20m The 
occurrence of an accompanying bright core-collapse 
supernova has bee n definitively ruled out for several 
short GRBs (e.g., iHiorth et ahl 120051; iKocevski et"all 
I20T0I) . 

Th e model of short GRBs as merge rs of comp act binary 
stars (jEichler et al.lll989tlPaczvhskil lT991: Nar avan et al.l 
II992D — either two neutron stars (NS-NS) or a neu- 
tron star and a black hole (NS-BH) — naturally pre- 
dicts old progenitors, and the observations discussed 
above appear to strongly support this model. However, 
other models permitting older populations do exist (e.g., 
iMacFadven et all |2005| ). and recent observations have 
painted a picture that appears more complex than it 
initially seemed after the discovery of the first few short- 
GRB hosts. A significant fraction of short GRBs with af- 
terglows actually have no clearly identifiable host (jBergeil 
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|2010D . suggesting that they have been flung large dis- 
tances from their original hosts or occur in optically un- 
derluminous, distant galaxies in surprisingly large num- 
bers. Finally, with the except ions of GRBs 05 0509B, 
050724, and possibly 100117A (Fon g et al.ll20lH ). every 
short GRB host to date (at this stage, a fairly formidable 
sample of >20 objects) has had at least m odest star for - 
mation in relation to its stellar mass ^e.g.. lBergei]l2010D . 
making it difficult to pin down a minimum progenitor age 
for the large majority of these events. The two "smoking 
gun" events of 2005 no longer appear to be particularly 
representative of the entire class. 

The relationship between short and long GRBs is 
at least vaguely analogous to that between Type la 
and Type II supernovae (e.g., IWheeleij Il981l) : Type 
la supernovae are often associated with moderate or 
old populations (and almost certainly are pro duced by 
compact objects, in this case white dwa rfs; iNomotol 
[TOSalBloom et al.|[20Tll: iNugent et al.ll201lD . while Type 
II supernovae are associated exclusively with young 
populations (and result from massive stellar core col- 
lapse). However, recent evidence suggests that Type 
la supernovae can actually be produced fairly rapidly 
in some cases, given their elevated rate in late- type 
hosts compared to what rnight be expected from an ex- 
clusively old population (iScannapieco fc Bildstcn 20()5i ; 
Mannucci et al.ll2005HSullivan et al.ll2006l: iLi et al.ll2011 : 
Maoz et al.ll2011|) — that is, the delay-time distribution 
is likely to be quite broad, with short-lived stellar pro- 
genitors as well as long-lived ones. 

It is possible that the short GRB delay-time distribu- 
tion may be similarly complex. Confirming this trend — 
and in particular, determining whether there is any need 
for a very young component associated with the youngest 
(and presumably most massive) stars, which coul d indi- 
catc a physically distinct pr ogenitor (e.g., Mctzgcr et al.l 
[20081: ILazzati et al.ll2010l: iVirgili et al.ll201lD — is of clear 
interest for better understanding the origins of the short- 
GRB population. 

Given the low short-GRB event rate (< 10 events 
yr"-*^) and uncertain selection biases affecting the pop- 
ulation (~ 25% of all Swift short GRBs are not de- 
tected by the XRT and cannot be localized), teasing 
such a tendency out of the full population is possible 
but challenging. Nevertheless, two recent studies have 
attempted to accom plish this using different techniques. 
IVirgili et all (|201lD analyzed the short-GRB redshift, flu- 
ence, and luminosity distributions compared to predic- 
tions from a delay-time convolved star- formation history 
of the Universe and concluded that a significant fraction 
of short GRBs may in fact originate from young stars. 
iLeibler fc BergeJ (|2010D used the properties of the host 
galaxies of a sample of short GRBs to estimate the age 
of the predominant stellar population, and similarly con- 
cluded that both short-delay and long-delay components 
are necessitated by the data. However, even the "short" 
component is consistent with a timescale of order 200 
Myr — much older than the lifetimes of massive stars 
that produce core-collapse supernovae, long GRBs, and 
nebular signatures in star-forming galaxies. 

The most unambiguous indicator of a very short-delay 
(<100 Myr) component would be the discovery of an- 
other "smoking gun" system, but this time with an 
extremely young stellar age. In particular, a short 



GRB within a starbursting galaxy whose current star- 
formation rate (SFR) is very large in comparison to its 
stellar mass (high specific SFR) would be a strong in- 
dicator that the event came from a star produced in 
the starburst, events whose characteristic times rarely 
exceed a few hundred Myr and are often much less 
(|Di Matteo et al.ll2008t [McQuinn et 111120091 12010I ). 

In this paper, we discuss the case of GRB 100206A 
as the first example of a short-duration burst whose 
host characteristics, at least on the surface, evoke such 
a system. In 5j2] we present our observations of the 
burst and its afterglow, showing it to be an unambigu- 
ous short-duration, hard-spectrum GRB with a faint, 
rapidly-fading afterglow detected only in the X-ray band 
despite deep, early optical imaging. We also present ob- 
servations of two galaxies in or near the X-ray error circle 
at a variety of wavelengths from optical through the mid- 
infrared. In Sj3] we analyze these data in further detail, 
measuring the metallicity and constraining the under- 
lying stellar population of the brighter host candidate, 
which we show to be a luminous infrared galaxy (LIRG) 
at redshift z — 0.4068 with a current SFR exceeding 
that of any previous short GRB by almost an order of 
magnitude. In 52] we argue that a posteriori statistical 
arguments strongly tie this galaxy to the short GRB, and 
examine the implications for the progenitor of the short- 
duration burst. While the current SFR is high, we note 
that the stellar mass is substantial and dominated by 
older stars, with no unambiguous evidence for starburst- 
ing, merging, or other short-lived features — consistent 
with an older progenitor, like that of other short-hard 
bursts. We summarize our conclusions in Sj5] 

2. OBSERVATIONS 

2.1. Swift-BAT and Fermi 

GRB 100206A triggered th e Burst Alert Telescope 
BAT: iBarthelmv et all |2005[) on the Swift satellite 
Gehrels et al.l 12004") at 13:30:05 on 2011 February 06 
(UT dates are used throughout this paper). BAT data 
are automatically red uced by our automated pipeline 
using the methods of iButler et~all (|2007[ ). The BAT 
light curve (Figure [Ij shows only a single spike start- 
ing at the trigger time and ending by 0.2 s, with no evi- 
dence of extended em ission (with Tgo = 0.200 ± 0.017 s; 
ISakamoto et al.ll2010t) . 

The GRB al so triggered t he Gamma-ray Burst Monitor 
(GBM; Meega n et aLll2009D onboard the Fermi satellite. 
The Fermi light curve is similar to that seen by Swift — a 
single bright spike lasting < 0.2 s (Tgo — 0.13 ± 0.05 s; 
lKienlinl l 201C)|). The spe c trum of the burst is quite hard; 
the best-fit " iBand et al.l ()1993f ) model indicates Epcak — 
^i9tll keV. 

These properties place GRB 100206A unambiguously 
within the short-durati on phenomenological c lass (see, 
for instance. Figure 1 of lLevesaue et al.ll2010cl) . 

2.2. Swift-XRT and UVOT 

Swift slewed immediately to the source and began 
pointed observation s with the X-Ray Telescope (XRT; 
iBurrows et"aIll2005D at 74.6 s after the BAT trigger, fol- 
lowed by obser vations with the Ultr aviolet Optical Tele- 
scope (UVOT; iRoming et all I2OO50 beginning at 78 s. 
Only a faint and rapidly fading X-ray afterglow is seen. 
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Fig. 1.— BAT light curve of GRB 100206A, taken from the sum 
of all four channel s of data (15-350 keV) using the methods of of 
IButler et a"Lll2007l . The main plot shows the light curve binned 
at 2 s (dark gray) and 0.2 s (light gray); the inset shows the light 
curve binned at 0.02 s. Times are referenced to t =949498223.86 sec 
(GPS). The burst is clearly short, with all detectable gamma-ray 
emission contained within a 0.25 s interval. There is no evidence 
of extended emission during the minutes after the trigger. 

dropping from 0.2 to 0.01 count s^^ during the first or- 
bit; in total, only 24 X-ray photons were detected by the 
instrument during this observation. Swift made several 
return visits, none of which resulted in a detection above 
the background level. 

During the first orbit, the event data are consistent 
with a light curve following a sharply falling, unbroken 
power law with decay index a = I.OIIq ^^. Given the 
limited number of counts and large Galactic hydrogen 
column (iVn ~ 10^^ cm~^), no clear statements can be 
made about the intrinsic X-ray spectrum. The best avail- 
able localization i s the UVOT-enhanced XRT position 
(|Goad et al.ll200l : a = 03h08"^39^03,(5 = -Kl3°09'25.3" 
(J2000, 3.3" uncertainty at 90% confidence). 

No detection is reported in UVOT observations at any 
epoch . UVOT upper limits are given by Krim m et al.l 
(|2010l ). In addition, we stacked all available UVOT it- 
band imaging of this object to try to provide a deep 
limit on any host-galaxy emission; the integrated limit 
is u > 21.8 mag. Unfortunately, given the relatively 
large Galactic extinction in the field {Eb-v — 0.382 
mag: [Schlegel et"aLlll998D . this value is not particularly 
constraining. 

2.3. Ground-Based Follow-Up Observations 

GRB 100206A was observed by many different ground- 
based i nstruments within the fir s t 24 hr aft e r the 
trigger dB hattacharva et al.' '2010'; 'Guziv et al.' '20101; 
Noda et a l. 2010; Lcvan et al. 2010; Loloudas ct al. 2010; 
Morgan et al.l 



Kuroda ct al. 2010; 

2010,; Mao et al. 201(3; Berger et al.ll^ _ _ , 

lAndreev et al.l l2010t iRumvantsev et al ■I l201flD . 



ft lYurkov et al. 

20ld iNicuesa et al.l 



None of these observations report a detection of a vary- 
ing source, although two stati c sources in t he XR T er- 
ror circle, i nitially noted by iMiller et al.l (|201Cl( l and 
iLevan et al.l (j2010. ). will be discussed below. Given the 
faint X-ray afterglow, this nondetection of variability is 
unsurprising, as explained in further detail in 

2.4. DeepSky 
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Fig. 2. — Unresolved images of the putative host galaxy (Gl) 
of GRB 100206A in the optical and near-infrared (NIR) from 
the 1.2 m Palomar Oschin Telescope (from the Palomar/DeepSky 
project), the 1.3 m PAIRITEL, and the WISE all-sky mission. Al- 
though these telescopes are relatively insensitive to typical galaxies 
at cosmological distances, a bright source centered just outside the 
XRT error circle is well detected in every filter except W4 (22 fim). 

We co-added 78 archival images from the DeepSky 
project at Palomar Observatory ([Nugent et al.l 120091) 
covering the field of GRB 100206A. The images were ob- 
tained in 2004-2008 from the Palomar-Quest Consortium 
at the Oschin Schmidt telescope. The limiting magnitude 
of the stack is i? w 23 mag. In the combined image we de- 
tect a faint (i? — 21.7±0.3 mag relative to nearby USNO 
catalog stars) extended source, centered slightly outside 
the current XRT error circle to its northeast (Figure [2]). 

The random appearance of a source this bright close 
to a small X-ray error circle is small, but not necessarily 
convincingly so (Pchancc ~ 0.05; see ^4.ip . However, this 
source has very unusual properties at other wavelengths 
that show it is far from being a typical galaxy. 

2.5. PAIRITEL 

The robotic Peter s Automatic Infrar ed Imaging Tele- 
scope (PAIRITEL; iBloom et all l2006f ) began observa- 
tions of the GRB 100206A field starting 13.1 hr after 
the burst in the J, H, and Kg filters simultaneously. No 
transient source was detected inside the error circle. 

However, the host-galaxy candidate first seen in the 
DeepSky archival images was well detected in all three 
filters. PAIRITEL detection of galaxies at cosmological 
distances f §2.9p is unusual, leadin g us to first speculat e 
that this galaxy may be a LIRG (jMorgan 
To improve the photometry of this source and verify the 
lack of fading behavior, we acquired additional imaging 
on 2011 October 20 and 2011 October 22, and combined 
data from all three epochs (total effective exposure time 
of 4.34 hr) with the exception of about 2 hr of if^-band 
data which was not usable due to poor sky conditions. 

Aperture photometry was performed using custom 
Python software, utilizin g Source Extractor (SExtractor; 
iBertin fc ArnoutsI 1 199(31 ) as a back end. Calibration was 
performed by determining the zcropoint for each image 
by comparison to 2MASS (Skrutskic ct al. 2006) magni- 
tudes using 19 calibration stars. As the galaxy is clearly 
extended in higher-resolution images, we employ a large 
aperture radius of 3". 

Indeed, no evidence of fading (to limits of J > 19.8, 
H > 19.3, Kg > 16.8 mag) is observed between the 
epochs (the Kg-hand constraint is poor due to the low 
quality of the second-epoch image). Photometry of the 
host galaxy using the final, combined stacks is presented 
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in Table m 

2.6. Gemini Imaging 

Two epochs of deep imaging were acquired with 
the Gemini Multi-Object Spectrometer (GMOS-N) on 
Gemini-North. In the first epoch (starting at 05:34, 14 hr 
after the trigger), five dithered exposures of 240 s each 
were obtained in the z band, immediately followed by a 
similar sequence of 5 x 240 s exposures in the i band. 
The i-band sequence was repeated five days later, be- 
tween 06:41 and 07:05 on 2010 February 12. Conditions 
were excellent during both observations, with an average 
seeing of about 0.8". 

The bright archival galaxy first seen in our DeepSky 
image is now resolved into an extended disk of about 5" 
diameter (Figure [3]). We denote this galaxy "Gl" in the 
remainder of the paper. A second, much fainter point- 
like source (which we denote "G2") is also seen ~ 8" to 
t he south of th i s obje ct . 

iBerger et al.l ()2010[ ). who triggered these observations, 
performed image subtraction between the two epochs of 
i-band imaging and reported no variation either between 
these epochs or in comparison with a WHT image taken 
7 hr after the trigger — including, in particular, no varia- 
tion of the po int-like sourc e G2 w hich had been initially 
suggested by iLevan et al.l ()2010[ ) as a candidate after- 
glow. (Indeed, our spectroscopy verifies that this object 
is a background galaxy; §2.91 ) 

We downloaded all Gemini frames from the Gemini 
Science archive and reduced them independently using 
the GMOS IRAF reduction tools. We first coadded the 
data obtained during each of the two separate i-band 
epochs, and subtracted the two images using HotPants^, 
confirming the lack of variability reported by Berger et 
al.^ We measure an improved 5(t upper limit of z > 25.9 
mag for any point source varying between the two frames 
over an 0.8" aperture within the XRT error circle. 

Given the lack of variation and comparable, good con- 
ditions over the two nights, we then stacked all ten i- 
band exposures over both epochs to produce a single, 
deep image. The five z-band frames from the first night 
were also coadded into a separate ^-band stack. Using 
secondary standard stars from our P60 calibration of the 
field ( i)2.8p . we performed aperture photometry of Gl 
using a radius of 3". Photometry of G2 was calculated 
using the same standards, using a 1" radius aperture; see 
Table □ 

2.7. Keck Imaging 

To obtain additional color information, on the night of 
2010-08-02 wc imaged the field of GRB 100206 A with the 
Low Resolution Imaging Spectrograph fLRIS: lOke et al.l 
Il995f ) on the 10 m Keck I telescope. Using the D560 
dichroic, we acquired images in the g and R filters simul- 
taneously, with 3 X 200 s integration in g and 3 x 180 s 
integration in R. Reduction was accomplished using a 
custom LRIS imaging pipeline. As with the Gemini im- 



°'http: //www. astro .Washington. edu/users/becker/hotpants .html . 

® One marginal source does appear with a nominal significance oT 
about 5tT at a = 03'^08™38.846=, & = -|-13°09'24.25" {J2000) but it 
seems to be mostly due to a weak negative artifact in the late-time 
reference image rather than the detection of positive variation in 
the first image. 




Fig. 3.— Deep imaging of the field of GRB 100206A from the 
Keck 10 m telescope and Gemini-North 8 m telescope, combined 
into a false-color image using the ij, i, and z filters. The bright 
galaxy Gl resolves into a highly inclined disk, with most of the 
southern half of the galaxy enclosed by the XRT error circle (shown 
in blue). A second, much more compact galaxy is also evident, 
centered at the southern edge of the XRT error circle, which we 
label G2. Isophotes of Gl show a subtle asymmetry that probably 
results from additional substructure within the galaxy (see also 
Figure [S]) The dashed green line indicates the position of the 
LRIS slit during our long-slit spectroscopy. North is up and East 
is left. 

ages, we performed aperture photometry of the two pos- 
sible host galaxies Gl and G2 using 3" and 1" radius 
apertures, respectively, using our secondary field stan- 
dards. Again, the photometry is reported in Table [T] 

2.8. P60 and Nickel Photometric Calibrations 

To calibrate the optical photometry, we acquired inde- 
pendent calibrations using the robot icized 60-inch tele- 
scope at Palomar Observatory fP60: [Cenko et al]|2006D 
and the Nickel 1 m reflector at Lick Observatory, both on 
2011 A ugust 09. Numerous standard fields from LandoI3 
(|2009D were observed throughout the night in the gRiz 
fihers (P60) or BVRI fihers (Nickel). The Nickel pho- 
tometry was then t ransformed to giz u sing the transfor- 
mation equations of lJester et al.l ()2005| ) . The two calibra- 
tions show good consistency, with < 0.04 mag systematic 
differences in each of the Riz filters for bright stars. A 
slightly larger offset of 0.1 mag is observed in the g band. 
To calibrate the galaxy photometry, we take the average 
of the two g, R, and i-band calibrations, adding a small 
calibration component to the uncertainty. Only the P60 
data are used to calibrate the z band. Final, calibrated 
photometry of the two galaxies near the XRT error circle 
is presented in Table [T] 

2.9. Keck Spectroscopy 

We acquired spectra using LRIS on the Keck I 10 m 
Jtelescope on two occasions, both using the D560 dichroic, 
the 600/4000 grism, the 400/8500 grating, and a 1" 
slit. The first epoch was obtained on 2010 February 07, 
shortly after the burst. Two exposures, each of dura- 
tion 600 s (on the red side; slightly longer on the blue 



The LIRG Host of Short GRB 100206A 



5 



TABLE 1 

Photometry of Galaxies Inside the XRT Error Circle of GRB 100206 A 



instrument 


Exp. time 


Filter 


Magnitude^ 


AB mag.'' 


Magnitude'' 


AB mag.'' 




(s) 




Gl (z = 


0.41) 


G2 (z = 


0.80) 


Keck I / LRIS 


600 


a 


23.83±0.17 


22.45 


26.55±0.40 


25.17 


Keck I / LRIS 


540 


R 


21.31±0.09 


20.51 


25.14±0.19 


24.34 


Gemini-N / GMOS-N 


2400 


i 


20.82±0.08 


20.06 


24.69±0.05 


23.93 


Gemini-N / GMOS-N 


1200 


z 


20.20±0.05 


19.63 


24.15±0.12 


23.57 


PAIRITEL 


16333 


J 


18.51±0.12 


19.06 


>19.7 


20.3 


PAIRITEL 


16333 


H 


17.25±0.09 


18.41 


>18.8 


19.9 


PAIRITEL 


8143 


Ks 


16.33±0.11 


18.03 


>17.7 


19.4 


WISE 




Wl 


15.74±0.06 


18.42 






WISE 




W2 


15.14±0.11 


18.47 






WISE 




W3 


11.23±0.16 


16.38 






WISE 




W4 


>8.58 


>15.18 







Observed value, not corrected for Galactic extinction. 
^ Corrected for Galactic extinction {Eb-v = 0.38 mag) 

side), were acquired at a slit position angle (PA) of 18.2° 
between 06:17 and 06:39. (LRIS is equipped with an 
Atmospheric Dispersion Corrector [Phillips et al. 200i|, 
enabling us to observe away from the parallactic angle ^ 
without slit losses.) Seeing conditions were poor 1.5") ^ 
throughout the integration. Nevertheless, the galaxy was 5. 
well detected in the red part of the spectrum, showing o 
strong emission lines of H a and [N II] , as we previously I 
noted (jCenko et al.ll2010l ). The slit orientation was for- | 
tuitously close to the galaxy major axis and these lines | 
show clear rotational structure. °* 

A second epoch was acquired on 2011 August 02, 
shortly following our imaging that same night ( ^2.7[ . A 
slit PA of 14.5° was used to simultaneously cover both 
Gl and G2 at the same time; one exposure was acquired 
of 900 s duration (on the red side; the integration on the 
blue side lasted an additional 10 s). This PA is still close 
to the major axis of the galaxy (Figure [3]). Conditions 
were excellent (0.8" seeing) and the spectrum is of signif- 
icantly higher quality than the one obtained during the 
earlier epoch, showing both lines of the [N II] doublet 
as well as [S II], H/3, and [O II]. Accordingly, only this 
spectrum is used in our final analysis. 

The spectra were reduced using custom reduction soft- 
ware written in IDL. The traces of the two galaxies 
were extracted separately using our custom software 
and flux-calibrated relat i ve to spectroscopic standards 
BD-F284211 (blue ; [Oki [T990l) and BD-l-174708 (red; 
lOke fc Gunnl [l98f). and by comparison of synthetic pho- 
tometry of the spectrum to the photometry from our im- 
ages. Both objects are observed to have multiple strong 
emission lines, identifying them as galaxies at redshifts 
of z = 0.4068 (Gl) and z = 0.803 (G2). In the case 
of Gl, strong rotational "shearing" is seen in the two- 
dimensional spectrum {v/c — 0.0014 between opposite 
ends of the galaxy, corresponding to 13 A in the vicin- 
ity of Ha), so a simple extraction would produce much 
poorer resolution than what is provided by the instru- 
ment. To remove this shear, we calculate the rotation 
curve of the galaxy by finding the maximum of the cross- 
correlation function along each row of the trace relative 
to the galaxy center in the region of the strong lines of the 
Ha-[N II] complex, and then resample each line (across 
the full trace) by the appropriate amount to remove the 
systematic rotation. This procedure allows us to achieve 
a resolution comparable to the instrumental resolution 
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Fig. 4. — Rotation curve for Gl, determined from Ha and [N II] 
emission lines. Note that the slit is 1" (5.4 kpc) wide and at a 
(slightly) different orientation than the major axis of the galaxy, 
so the true major-axis rotation velocities of the galaxy are likely 
to be somewhat larger than the values indicated here. (The galaxy 
itself is also likely inclined ~ 20° from edge-on; Figure [8]) The 
relevant region of the two-dimensional spectrum is shown in the 
inset. 

of 7 A. Following this step, the spectrum is extracted 
normally using a window of 3.24" (diameter). (The emis- 
sion lines of G2 are not resolved, and this spectrum is 
extracted normally with no additional steps.) The final 
ID spectra of Gl and G2 are presented in Figures [5] and 
m respectively. 

2.10. WISE Archival Observations 

We se arched the Wide-Fi eld Infrared Survey Explorer 
(WISE; IWright et all l2010l ) archive for observations at 
the position of this GRB. The field was covered by the 
WISE preliminary data release, and a source is clearly 
detected at the position of Gl in both the archival im- 
ages and in the photometric catalog in filters W1~W3 
(no detection is evident in W4), as shown in Figure [5J 
In Table [1] we present catalog magnitudes of this object 
from the archive. 

3. ANALYSIS 

3.1. Afterglow and Kilonova Models 

Despite extensive optical and NIR follow-up observa- 
tions during the first night (including some data on time 
scales as short as minutes and a very deep image with 
Gemini), no afterglow was detected from this event. Us- 
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Fig. 5. — Spectroscopy of Gl, the putative host galaxy. The combined LRIS blue+red spectrum is shown in the main panel, with 
the positions of detected emission lines indicated. In addition, the red curve shows the best stellar-continuum model from a fit to our 
broadband (and synthetic narrowband) photometry. The small deviation between this model and the observations at >9500 A is probably 
due to uncertainties in the spectrophotometry in the long-wavelength region. At top, the insets show regions around specific lines, including 
[O II], the Balmer absorption fines, H^-|-[0 III], and Hq:-|-[N II] + [S II]. 



2000 



3000 



Rest Wavelength (A) 
4000 



5000 




4000 



5000 



6000 



7000^ 
Wavelength (A) 



8000 



9000 



10000 



Fig. 6.- 
show this 



— Spectrum 
to be a star- 



of G2, 
formint 



the faint galaxy at the south end of of the XRT error circle. Secure detection of the [O II] and [O III] lines 
5 galaxy at z = 0.803. The vertical gray lines indicate the centers of strong sky emission lines. 



ing standard assumptions about the intrinsic spectral in- 
dex of a GRB afterglow, we examined whether these Um- 
its might usefully constrain the properties of the GRB or 
the extinction column. 

In Figure [7] we plot the "light curve" of the X-ray 
afterglow (consisting of two binned detections followed 
by upper limits) and scale all reported UV, optical, or 
NIR upper limits from the GCN Circulars to the X-ray 
band (1 keV equivalent photon energy) based on an as- 
sumed afterglow spectral index of /3ox = 0.5, the mini- 
mum value expected in the synchrotron afterglow model 
(jJakobsson et al.l 120041 : a value below this line would in- 
dicate a "dark" burst). Any upper limits lying below the 
X-ray light curve are inconsistent with this basic model 
and would require either an unusually blue/hard intrin- 
sic spectral index or extinction within the host galaxy to 
suppress the optical flux. As is evident from the figure, 
none of the limits constrain the afterglow in this way. In 
fact, given that the second X-ray point is a nondetection, 
even if we assume a much softer/redder spectral index 



(up to f3ox = 1-1, which corresponds to an unbroken 
spectral index between the X-ray and optical given typ- 
ical X-ray spectral indices and represents the maximum 
value expected in the synchrotron model), the optical 
photometry still imposes no constraint on any additional 
host extinction. 

The time of the deep Gemin i limit is better timed 
to constrain emission from a iLi fc Pac zvhski' (1998|) 
mini-supernova or "kilonova" (see also iKulkarni .2005t 
iMetzger et al.|[2010l ). which is expected to peak at t « 1 
day. After correcting for Galactic extinction, the GIVIOS 
i-band limit corresponds to a specific luminosity limit of 
vL^ < 10*2 erg s"i at i = 11 hr, if the GRB is at the red- 
shift of the brighter galaxy (z = 0.4068) and is unextin- 
guished by dust with in its host. This i s com parable to the 
limits presented by iKocevski et al.l (|2010t ) and broadly 
rules out models with large energy conversion factors and 
ejecta masses (/ > 10"^ and M > IQ-'^-^Mq), but is 
significantly less constraining than the limits presented 
on radioactive-powered emission from GRB 050509B by 
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Fig. 7. — X-ray and UV/optical/NIR observations of GRB 
100206A. The X-ray light curve is plotted as flux density {fi,) at 
1 keV. To meaningfully place the UV/optical/NIR observations (all 
of which are upper limits) on the same plot, we calculate the flux 
corrected for Galactic extinction and extrapolate it into the X-ray 
band assuming /3 = 0.5 {f^, = u~^). Given the extremely faint 
X-ray afterglow, none of the optical limits is constraining. 



iHiorth et al] H2005f ). Furthermore, even our weak hmit 
is subject to the strong caveat that the extinction within 
this host galaxy appears to be quite large f i^S.Sp and the 
assumption of no extinction may not be valid. 

3.2. Morphology 

The galaxy Gl has a disk-like morphology (Figure [3]), 
and is seen at a high inclination. It is also quite large, 
with significant emission extending out 1.5" (8 kpc) in 
each direction away from the nucleus. Some subtle asym- 
metry is visible: the location of maximum flux is dis- 
placed by about one pixel (0.15") from the center of the 
outer isophotes, and there is a hint of displacement be- 
tween the northern and southern sides of the disk from 
the nucleus, possibly suggestive of spiral or bar structure. 

To gain further insight into the struct ure of the galaxy , 
we used the software package GALFI T (iPeng et al ll2002D 
to model the system as an inclined ISersid ()1963() disk. 



Our modeling is performed using the Gemini GMOS i- 
band image, which has the best signal-to-noise ratio and 
seeing quality of the filters used. A bright, isolated, non- 
saturated, nearby star in the image is used as the point- 
spread function (PSF) model. Because a simple Sersic 
fit gives significant residuals (as expected, given the vi- 
sual asymmetry), we extend this basic model by succes- 
sively adding individual Gaussians initially centered at 
locations of positive residuals and iteratively repeat the 
fit until no residuals are evident above the background 
noise. We found three such components to be needed 
within the physical extent of the disk: one each on either 
side of the nucleus in the disk plane, each at a displace- 
ment of 1.3" (7 kpc in projection), plus an additional 
component displaced 1.1" east (6 kpc) of the nucleus. 
(A fourth component 3.2" northeast of the galaxy is also 
required by the fitting procedure, but is probably not 
related to the system.) These additional components 
combined only contribute ~ 20% of the total fiux of the 
system, which is still dominated by the Sersic disk. 
Whether or not these additional parameters are added. 




^ 5.4 kpc 




Gemini-N/GMOS i-band 



Deconvolved 



Fig. 8.— The field of GRB 100206A imaged in the i band us- 
ing the Gemini-North 8 m telescope, with no additional processing 
(left) and as deconvolved by our GALFIT model (right). In addi- 
tion to a Sersic exponential disk, three other compact components 
are required to obtain a fit without residuals: two along the major 
axis displaced 7 kpc in projection from the center and an additional 
object close to the minor axis. 

similar results for the Sersic parameters are achieved 
(Sersic index n — 0.58, and axis ratio b/a = 0.31 for a fit 
with the disk only and no extra components; n = 1.27, 
b/a = 0.40 for the final model with four extra compo- 
nents). Evidently, the stellar light in this system follows 
a nearly exponential disk (n = 1; typical for disk galax- 
ies), inclined at an angle of ~ 70° from face-on. 

The decomposition described above also functions as 
a basic deconvolution of the image near the galaxy. In 
Figure [8] we show both the observed i-band image and 
its deconvolved equivalent, which is the sum of the model 
produced by GALFIT {not re-convolved with the PSF) 
and the fit residuals. The residuals are added both to 
accurately show the noise level, and also to avoid falsely 
removing any real signal that may be left in the residuals. 

Deconvolution should always be approached with cau- 
tion, and our modeling cannot identify faint, small-scale 
(^ 1") structures that are likely to be present or more 
complex structures such as bars or arcs. Nevertheless, 
since the observed components are all separated on scales 
larger than the 0.8" seeing disk, we expect that they 
likely represent real substructure in the galaxy — pro- 
jected spiral arms or localized intense star-forming re- 
gions (the regular Ha rotation curve in Figure |4] sug- 
gests they are not merging galaxies). Further observa- 
tions with better image quality (such as from the Hubble 
Space Telescope) would be necessary to unambiguously 
resolve the structure of this system. Nevertheless, it is 
clear that the host is a morphologically complex but pre- 
dominately disk-like galaxy. 

3.3. Emission-Line Measurements 

The clear detections of numerous, bright emission lines 
in the spectrum of Gl allow us to apply standard diag- 
nostics of the SFR, extinction, and metallicity to this 
galaxy. 

We first extract the emission-line fluxes by assuming 
a flat continuum with no underlying stellar absorption 
component. Taking a small region around each emission 
line or line complex, we fit a Gaussian function (or, in the 
case of blended lines, several summed Gaussians) added 
to a linear component representing the underlying flux 
with the mpf it package within IDL. For each line, the 
integrated fiux and uncertainty are obtained from the fit. 
The results are presented in the third column of Table [2j 

While the strong Ha line (combined with the mod- 
est [N II] /Ha and [O IIIJ/H/3 ratios which associate 
this line with star formation rather than AGN activ- 
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TABLE 2 

Fluxes of Detected Emission Lines in G1 and G2 



Galaxy Line Species 


Jrtest Waveiengtii 


Unsubtractea riux 


oUDiracT^ea r lux 






(A) 


(10 erg s ) 


(10 erg s ) 


Gl 


O II]= 


3727 


16.2 ± 1.8 


19.6 ±2.1 


H/3 


4861.33 


12.7 ± 1.8 


21.8 ± 1.9 


[O III] 


5006.84 


5.9 ±1.5 


5.1 ± 1.2 


Ha 


6562.82 


108.6 ± 1.4 


112.8 ± 1.6 




N II] 


6548.06 


12.8 ± 1.1 


13.1 ± 1.2 




N II] 


6583.57 


46.4 ± 1.2 


45.7 ± 1.3 




S II] 


6716.44 


20.4 ±3.2 


20.4 ±3.1 




S II] 


6730.82 


14.1 ± 1.8 


13.6 ± 1.7 


G2 


O II|= 


3727 


6.0 ±0.8 




H/3 


4861.33 


2.8 ±1.7 




[O III] 


5006.84 


3.1 ±0.4 





^ Assuming a flat continuum, corrected for Galactic extinction. 
^ After subtracting the model continuum, corrected for Galactic extinction. 
Combined flux of both lines of the [O II] doublet. 



ity; [Baldwin et al] |1981[ ) gives unambiguous evidence 
of rapid, recent star formation, several other indicators 
suggest that an older stellar component is also present. 
Specifically, while many of the higher-order Balmer lines 
fall in the wavelength region of the dichroic where sensi- 
tivity is relatively poor (~ 5550-5790 A), we do see a dip 
at the location of H(5 and lower-significance dips at the 
locations of H?7 and H^, as well as (possibly) Ca H&K 
(inset of Figure [5]) . Interestingly, we see neither emission 
nor absorption at the location of H7 where sensitivity is 
good, suggesting that the emission and absorption are of 
comparable strength, resulting in a fiat apparent spec- 
trum (the instrument resolution of 7 A is comparable 
to the intrinsic width of the Balmer line). We also see 
evidence for a strong Balmer break in the photometry 

(SSD- 

The presence of Balmer absorption features raises the 
prospect that the simple procedure above used to mea- 
sure the emission-line fiuxes may be systematically un- 
derestimating some line strengths. The large equivalent 
width of Ha suggests that it is unlikely to be affected by 
uncertainties in the underlying continuum by more than 
a few percent, but the weaker H/3 line may be much more 
strongly affected, which would impact the spectroscop- 
ically derived extinction (and therefore the extinction- 
corrected SFR). Consequently, after deriving an estimate 
of the stellar continuum from population synthesis mod- 
eling ( ij3.51 below), we subtract the continuum fiux given 
by that model and then repeat the line fits with the con- 
tinuum removed. The results are given in the fourth 
column of Table [2J 

For galaxy G2, we only present the line fiuxes using 
the basic linear continuum-subtraction model. As the 
continuum trace of this galaxy is extremely weak, this 
is a good assumption at the level of uncertainty in our 
measurements. The fiuxes of the probable lines (note 
that H/3 is blended with a strong night-sky line and not 
clearly present) are given in Table [21 

3.4. Chemical Abundance Analysis 

The strong and unambiguous detection of numerous 
emission lines makes it possible to measure the metal- 
licity of this short-GRB host galaxy; this has been suc- 
cessfully done for only a few other short-GRB hosts so 
far. IVIetallicity is expected to infiuence not only the 



lives of stars but also the outcome of their deaths as dif- 
ferent kinds of explosion; thus, the burgeoning field of 
metallicity studies for different kinds of expl osive events 
from both observati onal (for a review se e iModiadl201lD 
and t heoretical fe.g.. iHirschi et alll2005t lYoon fc Langeil 
120051) perspectives has emerged over the last few years. 

The nebular oxygen abundance 12 + log(0/H) is the 
canonical choice of metallicity indicator for studies of 
the interstellar medium. Using the derived line fiuxes 
from Table [5] and correcting for host-galaxy reddening 
Eb-v = 0.60 ± 0.42 mag assuming Ha/H/3 = 2.86 from 
Case B recombination (Osterbrock 1 989t ). we measure the 
nebular oxygen abundance of the host ga l axy fo l lowin g 
the technique descri bed bv iModiaz et all |2008l 120111) . 
Using the scales from lPettini fc Pagell (|20o4 fPP04). we 
find 12 + log(0/H) = 8.74;° °2 based on the [N II]/Ha 
diagnostic (PP04-N2), and 12 + log(0/H) = 8.8ltg:i^ 
based on the [O III]/[N II](PP04-O3N2) prescription. 
Likew ise, we find 12 + log( 0/H) = 9.17to °i on the scale 
of Kewlev fc Dopital ([20 ^ (KD02) and 1 2 + log(0/H) 
= 8.82;°; on the scale of llMcGaughl (Il99l (1V191). Here 
we have computed the uncertainties in the measured 
metallicities by explicitly including the statistical uncer- 
tainties of the line-fiux measurements, a conservative es- 
timate of the uncertainty arising from continuum absorp- 
tion contamination, and those in the derived host-galaxy 
reddening, and we propagate them into the metallicity 
determination. The ionization parameter g, iteratively 
derived in the KD02 models, is logio(q) = 7.59j;Q'27. 

Given the most recent estimate of the solar oxygen 
abundance [12 + log(0/H) = 8.69: lAsplund et al.ll2009| . 
together with the r elative robustnes s of the Tg-based 
metaUicity scale (Bres olin et"alll2009[ ) . which the PP04- 
03N2 scale is close to, we conclude the true value of 
the host metallicity Z is most likely to be closer to 
the lower end of the estimated range, approximately 
12 + log(0/H) = 8.8, or Z = 1.1 Z©; mildly super- 
solar. However, most pr evious metal licity work on short 
GRB host galaxies (e.g.,'BergeiJl2009') has used the KD02 
scale exclusively, which tends to produce systematically 
higher values than PP04 (see detailed discu ssions in 
iKewlev fc Elhsonl l2008t iMoustakas et all l2010l and ref- 
erences for possible reasons for the systematic offsets 
between different abundance diagnostics). For more di- 
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rect comparison with other GRB host galaxies, then, the 
KD02 value of 12 + log(0/H) = 9.17j:j]:?? (a significantly 
larger value of Z = SZq) is more informative. 

In either case, we conclude that Gl is the most 
metal-rich host galaxy of any GRB (long or short) 
to date. To our knowledge, the only super-solar 
long GRB host galaxies (as determined via emission 
spectroscopy, via the KD02 scale^°) a re the hosts of 
GRB 20819 (12 -I- log(0/H) = 9.0±0.1: lLe^esaue et al.l 
20T?)bh. G RB 050826 (12 + log(0/H) = 8.83 ± 0.1; 
Levesaue et al..,2010ai) an d GRB 051022 (12 + log(0/H) 
8.77: iGraham et all l2010l 1201 1[ ). Among short 
GRBs, the host galaxies of GRBs 051221A, 061210, and 
070724 A aU have me tallicit i es of 12 + log(0/H) = 8 8-8.9 
(also using KD02; IBergeij l2009t ) iKocevski et all ()2010f ) 
report a KD02 value of 12 + log(0/H) = 9.1 for GRB 
070724A. Using the same scale, our measured value for 
GRB 100206A (12 + log(0/H) = 9.17^^?) identifies 
this galaxy as the most metal-rich GRB host known. 

3.5. Spectral Energy Distribution 

We have recently developed a simple, flexible code 
for fitting multiwavelength photometric observations of 
galaxies. Implemented in IDL using mpfit, the code 
uses a small number of smoothly varying fundamental 
parameters (metallicity, mass, current SFR, and param- 
eterizations of the past star-formatio n history) and a grid 
of pop ulation-synthesis models from iBruzual fc CharlotI 
(|2003f l to calculate the stellar spectral energy distribu- 
tion (SED). This is extinguished by dust using one of 
several standard extinction laws, and the absorbed en- 
ergy is reradiated in the form of a multi-temperature 
dust graybody and mid-IR polycyclic aromatic hydrocar- 
bon (PAH) emission features. Emission lines are also in- 
cluded, scaled using the relations from Kcnnicutt (1998), 
iKewlev fc Dopital |2002^, and lKewlev et al.l (|2004i) using 
the current SFR, metallicity, and ionization. 

We initially fit our models using only the broadband 
photometry, fixing the SFR and extinction using the 
Ha and H/3 emission-line flux measurements. However, 
the probability of strong, unresolved stellar absorption 
underlying H/3 makes it difficult to apply these con- 
straints in a self-consistent way — and, furthermore, 
the age/extinction degeneracy makes it nearly impossible 
to place useful constraints on the stellar population age 
without photometry in the vicinity of the Balmer and 
4000 A breaks. Fortunately, our flux-calibrated spec- 
trum is of sufficient quality to fill these gaps. Using this 
spectrum, we calculate synthetic photometry using a se- 
ries of customized "filters" covering several independent, 
critical regions of the spectrum — the major absorption 
and emission lines ([O II], [O HI] A5007, [N II] A6854, 
[S II] AA6717, 6731, H<5, H7, H/3, Ha) as weU as inter- 
line continuum regions just blueward of the Balmer break 
(5650-5730 A and 5850-6050 A). (With the exception of 
the [O II] line, we do not attempt this technique blue- 
ward of the Balmer break due to the weak continuum 
trace and relatively large uncertainty in the overall flux 
calibration of the blue side.) We exclude the W3 fllter 

Among the hosts discussed in this paragraph, only GRB 
0208 19 was also measured o n the PP04 scale: the reported value 
from lLevesaue et al.l l|2010bl ') is 8.7 ±0.1, similar to what we derive 
for GRB 100206A on the same scale. 



in these flts, since the model is subject to large system- 
atic uncertainties in this wavelength region regarding the 
fraction of dust emission in PAH lines and the possible 
presence of highly embedded star formation. (However, 
all of our models which reasonably flt the optical/NIR 
points also accurately predicted the W3 flux within 2cr, 
and including this fllter does not qualitatively change any 
results.) 

Three different models of the past star-formation his- 
tory were attempted (see inset of Figure IH]): purely con- 
stant (up to and including the present-day value), con- 
stant (with an instantaneous change ai t = 10'' yr), 
and exponentially falling with a 100 Myr e-folding time 
(also with a step at t = 10^ yr). For each of these 
star-formation histories, we attempted two models of 
the dust attenuation: either a single dust screen applied 
to all stars uniformly, or with two different screens — 
one applied to the young (< 100 Myr) stars and neb- 
ular lines, the other to olde r stars (> 100 My r). We 
try Si nall Magellanic Clou d (|Gordon et al .''2003V Milky 
WayJCardelh et al.llT989l) . and Calzetti (Calzct ti et al.l 
l2000r ) extinction "laws" for each fit (Local-Group type 
extinction as implemented using the gen eralized param- 
eterization of iFitzpatrick fc MassallT990[ ) and report the 
result producing the best in Table ID 

With the exception of the purely constant, single- 
extinction model (which produces a poor fit), all of these 
models produce reasonable fits to the combined broad- 
band and pseudo- narrowband photometry of Gl. These 
different models are shown in Figure [9l with the param- 
eters outlined in Table Different assumptions natu- 
rally lead to an intrinsic dispersion of properties, but es- 
sentially all good fits share several features in common. 
The current SFR is large but not extreme (20-40 M©), 
and the stellar mass is also quite high, (3-7) xlO^" M©. 
The (maximum) age of the stellar population is at least 
0.5 Gyr (at least 1 Gyr for the constant star-formation 
history). The extinction is also large. Ay = 1.8 ± 0.2 
mag (for a single extinction) or Ay^young = 2.7 ± 0.4 
mag and Ay.oid = 1-1 i 0.2 mag (age-dependent extinc- 
tion). All models require a large bolometric luminosity 
of L = (3 — 4) X IO^-'^Lq, the large majority of which is 
emitted in the far-IR, thus classifying the galaxy as an 
unambiguous LIRG. As we do not have any far-IR mea- 
surements directly sampling the thermal dust emission, 
it is possible that the luminosity is even higher than this 
if an additional deeply embedded {Ay > 50 mag) com- 
ponent is present in the galaxy. However, the relatively 
modest mid-IR emission in the W3 filter and the W4 
nondetection suggest that such an extra component is 
not likely to be present, as is the case for most high- 
redshift LIRGs (Reddy et al. 20 10|). but unlike nearbv 
LIRGs and ULIRGs fe.g.. Svmeo^dis et al.llMsl) . 

G2 is only detected in four filters, and no hydrogen 
lines are evident in the spectrum, so it is not possible 
to constrain the properties of the galaxy in any detail. 
Assuming no extinction, the galaxy is consistent with 
being a low-mass, non-starbursting galaxy with a modest 
current SFR of -0.5 Mq yr^^ and mass 4 x lO^M©. 

4. DISCUSSION 

4.1. Association of Gl with GRB 100206 A 

We have previously suggested that the probability of 
chance alignment of an IR-bright source with an XRT po- 



10 



Perley et al. 



TABLE 3 

Gl Metallicity Determinations 



Method Type'' Relevant lines'^ 12 + log(0/H) Z/Zq'^ 



PP04-N2 


Empirical 


[N II], Ha 


74+0.02 


1 19+0.05 


PP04-O3N2 


Empirical 


[O III], H^, [N II], Ha 


01+0.10 

°-°-'--0.12 


1 09+0.34 

^•■^^-0.32 


M91 (R23) 


Theoretical 


[O II], [O HI], ni3 


09+0.11 

°-°^-0.16 


1 qc; + 0.87 

1.00_Q 3Q 


KD02 (combined) 


Theoretical 


[O II], [O HI], [N II], [S II], H/3 


g 1 7+0.09 


(19+0.87 
■'•^^-0.93 


Fitted^ 


Theoretical 


[O II], [O HI], [N II], [S II], Ha, R/S 


9.04 ± 0.09 


2.25 ± 0.4 



Principle of the metallicity calibration: theoretical (using photo-ionizaion models) or em- 
pirical (calibrated to observations of local star- forming regions.) 



^ Lines directly used as part of the metallicity determination prescription. Most models 
additionally include Ha and H/3 to determine the host extinction, which is then used to correct 
the metal line fluxes. 

Gas-phase oxygen abundance (in Solar units, for log[0/H]Q = 8.69). 
Weighted average of the lower three models in Tabled see ^3.5l for details. 



TABLE 4 

Results of Model Fits to Photometry and Spectroscopy of Gl 





SFRo'' SFRav'^ 


Age"^ 




Mass' 


Dusts 


Av (old)'' 


Av (young)' 






(M0 yr-i) (M0 yr"!) 


(Gyr) 


(Myr) 


(109 M0) 


law 


(mag) 


(mag) 




2.56 ±0.51 


11.6 ± 1.8 


11.5 ±3.3 


00 


133 ± 31 


smc 


1.2 ±0.1 


25.47 / 14 


2.16 ±0.46 


17.0 ±3.0 41.0 ±39.4 


0.9 ±0.8 


CXD 


36 ± 14 


smc 


l.h 


±0.2 


17.00 / 13 


2.58 ±0.36 


19.4 ±4.2 94.3 ±26.6 


0.6 ±0.2 


100 


54 ±5 


calz 


l.t 


±0.2 


13.48 / 13 


2.32 ±0.36 


24.0 ±6.3 


2.4 ±0.9 


00 


57 ±16 


calz 


1.3 ±0.1 


2.2 ±0.3 


12.82 / 13 


1.94 ±0.75 


37.4 ±13.1 18.9 ±14.0 


3.4 ±2.3 


00 


65 ±21 


calz 


1.1±0.2 


2.7 ±0.4 


9.16 / 12 


2.27 ±0.46 


34.1 ± 14.6 82.0 ±24.4 


0.6 ±0.2 


100 


52 ±5 


calz 


1.6 ±0.3 


2.5 ±0.5 


11.88 / 12 


^ Gas-phase oxygen 


abundance (in Solar units, for lo^ 


;[0/H]o = S 


.69). 













^ Current SFR. 

Average SFR. 

Age of formation. 
^ SFR decay time scale (for exponential-decline model). 
^ Total stellar mass. 

^ Best-fit dust extinction curve (Milky Way, Small Magellanic Cloud, or Calzetti). 
^ Extinction for t > 100 Myr stars. 
^ Extinction for t < 100 Myr stars. 



TABLE 5 

Chance Alignment Parameters for 01 and G2 



Galaxy 


Filter 






-^chance 


-^chance, all SHBs 






(mag) 


(deg~^ mag^') 






Gl 


a 


22.33 


5730 


0.05 


0.80 




R 


20.38 


3160 


0.03 


0.58 




i 


20.09 


2280 


0.02 


0.46 




z 


19.67 


2280 


0.02 


0.47 




J 


18.16 


2000 


0.02 


0.42 




H 


17.03 


820 


0.007 


0.20 




Ks 


16.19 


740 


0.006 


0.18 




3.6 ^m 


15.74 


860 


0.008 


0.21 




4.5 ^m 


15.14 


950 


0.008 


0.23 




W3 


11.23 


170 


0.0015 


0.04 


G2 


a 


25.21 


1.1 X 10^ 


0.61 


1.0 




R 


24.19 


7.5 X 10'' 


0.48 


1.0 




i 


23.95 


8.3 X 10'' 


0.52 


1.0 




z 


23.62 


1.0 X 10^ 


0.58 


1.0 
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Rest Wavelength ((J,m) 




Wavelength (\im) 

Fig. 9. — Photometry of the putative host galaxy Gl of GRB 
100206A, fit with stellar po pulation m odels using our own imple- 
mentation of the Bruzual fc Charloti | |2003l ) templates (including 
nebular lines, dust extinction, and mid/far-IR dust emission) for 
different assumptions of the star formation history. The broadband 
photometry is supplemented by synthetic narrowband photometry 
of major emission and absorption-line regions (and interline re- 
gions near the Balmer break). Three different star- formation his- 
tory models are shown, all of which produce similar results (see 
also Table [4t . The green curve shows a strictly continuous star- 
formation history from the formation of the galaxy until the present 
time. The blue curve is also constant, except for an instantaneous 
change in the recent past. The red curve assumes an impulsive 
star-formation episode at some point in the past with exponential 
decay time t = 200 Myr. All three models use an age-dependent 
dust screen. Broadband photometry is indicated with large yellow 
points; synthetic narrowband photometry is indicated with smaller 
points. Empty colored squares show the synthetic fluxes for each 
filter. The gray lines show contours of constant AB magnitude. 

sition is very low. Here we quantify this calculation, and 
demonstrate that the case for associating GRB 100206A 
with the galaxy Gl is quite strong. 

Figure 1 shows the current XRT error circle relative to 
our optical imaging. Note that while the galaxy is cen- 
tered outside the XRT error circle, a significant amount 
of flux from the southern end of the disk is in fact con- 
tained within the X RT error circl e. 

Roughly following iBloom et alj ()2002|) , an estimate of 
the probability of chance association Pchance can be pro- 
vided by 

-^chance — 1 CXp( ^assoc ^ ^m<T7iobj)' 

Here, ^assoc is the area on the sky of the region in which 
an error circle centered in that region would still lead us 
to connect the GRB and the galaxy. We take this as a 
circle of radius equal to the sum of the galaxy's major 
axis and the 90% confidence radius of the XRT position: 
~ 6", so Aassoc = 113 square arcseconds. 

The term nm<mahi density on the sky of objects 

of equal or greater brightness. Because of Gl's unusual 
colors, this factor is strongly dependent on the filter 
chosen. Rather than restricting our calculation to a 
single filter, then, we calculate Pchanco for all filters in 
which the galaxy was detected. The number density 
of galaxies on the sky for various filters was t aken 
from a variety of sources :_ lYasuda et alJ (|2Q01[ ) for 
the SDSS _giz bands; H ogg et al l (I1997D for the R 

(|2006[ ): llmai et al.l 
the J, H, and K 
the number-count 



plots in iJarrett et al.l ()2011l ). while for the Wl and 
W2 filters we adopt the more p recisely det e rmine d 
Spitzer/IRAC number counts in iFazio et al.l ()2004[ ). 
transforming the WISE filters to equivalent IRAC fluxes 
using synthetic photometry of our model. Results were 
checked for consistency with the compilation graphs at 
http://astro.dur.ac.uk/~nm/pubhtml/counts/counts.html 



band; iJarrettI (l20b4l\: IFrith et al 



(|200l :" lMaihara et al.l (l200l for 
bands. For the W3 filter we use 



for relevant filters. 

The sky densities and Pchancc values are presented in 
Table [51 The fourth column lists the basic Pchancc as de- 
scribed above, the value of which is small almost regard- 
less of the filter examined — ranging from a few percent 
in the optical filters down to ~ 10"'^ in the W3 filter, 
similar to the strength of ass ociation of GRB 05 0509B 
to its putative elliptical host ("B loom et al.l 120061 which 
was based on r-band photometry). 

To adopt a more skeptical perspective, it should be 
noted that the number of localized short GRBs is no 
longer small, and curious chance alignments are bound 
to happen eventually. We therefore also calculate the 
Reliance using the combined area covered by all 46 short 
GRBs with X-ray localizations to date — ~ 733 square 
arcseconds, or (adding the 3" major-axis of the galaxy) 
3600 square arcseconds. The value of Pchance is 0.18- 
0.8 for the optical and NIR filters, an obviously less 
convincing result, indicating that the global short-GRB 
sample size is indeed approaching a point where bright 
foreground or background interloper galaxies can be ex- 
pected somewhere in the full sample. However, when 
comparing to the W3 filter magnitude — where the un- 
usual characteristics of this galaxy are most pronounced 
— the probability is still quite low at only 0.04. In other 
words, the current short-GRB sample size is still a factor 
of ~ 5-10 too small for chance alignment with a compa- 
rably mid-IR bright galaxy to be expected. 

An alternative hypothesis that GRB 100206A orig i- 
natcs from G2 was first presented bv lBerger et al.l ()2010[) . 
While this certainly cannot be ruled out, we judge it to 
be relatively unlikely given the faintness of this galaxy 
and the size of the XRT error circle: for all four filters, 
Reliance of the individual association is ~ 0.5 (i.e., a ran- 
domly positioned error circle of this size will enclose a 
comparable source approximately half the time). 

Based on these considerations, we argue that the host 
galaxy of this event is Gl. As with all short-GRB host 
associations to date, this one is based entirely on a pos- 
teriori arguments and — given the absence of an absorp- 
tion spectrum — is impossible to prove conclusively. If 
this GRB is associated, however, the unique properties 
of this galaxy (among short-GRB hosts so far) have the 
potential to impact our understanding of the short-GRB 
progenitor in important ways. In the remainder of the 
paper, we will assume that this association is correct and 
examine its implications. 

4.2. Characteristics of the Galaxy and Implications for 
the Progenitor 

With an SFR of at least 15 Mq yr^^ (and 30-40 Mg 
yr~^ in our best-fit models), Gl is the most rapidly star- 
forming host of a short-hard GRB to date. The next- 
highest reported SFR among short-GRB h osts in the lit - 
erature is only 6.1 M© yr"! (GRB 06080 l;|Berge3[2009^, 
much less than our inferred value. (After that, the next- 
highest short-GRB SFRs belong to GRBs 061217 and 
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070724, both at a mere 2.5 M© yr"^)" The intense on- 
going star formation raises the possibihty that the GRB 
itself may be the product of a short-hved star, coming 
from a young stellar progenitor produced at the same 
time as the large population of massive stars currently 
powering the bright nebular and mid-IR emission. How- 
ever, more detailed analysis of the galaxy suggests that a 
large population of evolved stars must be present as well, 
and that the current SFR is not even unusual relative to 
the rest of the galaxy's history. 

Unfortunately, there is no direct means of determining 
which of these two populations produced the GRB. The 
X-ray position does not locate the burst within the host 
(and, in fact, leaves open the possibility of a halo ori- 
gin), nor can any properties of the absorbing column be 
inferred from the afterglow. 

Bereft of any direct means of evaluating the age of 
the progenitor — and, thus far, only a single example of 
a dusty, rapidly star forming, short-GRB host to date — 
we must resort again to statistical arguments. Relative 
to the Universe as a whole, is this galaxy more likely to 
produce explosions from old progenitors, or young ones? 
To phrase this question more precisely (under the sim- 
plifying assumption that the short-GRB rate is directly 
proportional to the number of stars of the appropriate 
age range and ignoring secondary effects such as detec- 
tion biases), we can ask the following: as a relative frac- 
tion of all stars in the Universe at that redshift, are there 
more "old" stars or "young" stars in LIRG-like hosts at 
z « 0.4 ? (For now, we will encapsulate the "unusual" 
characteristics of this host galaxy entirely in its designa- 
tion as a LIRG.) Fortunately, this type of question has 
been addressed in great detail by recent studies: the sites 
of star formation, the build-up of stellar mass, and the 
role of the LIRG phase in galaxy evolution have been 
major foci of recent work in extragalactic astronomy. 

We will first address the fraction of star formation oc- 
curring in LIRGs. At higher redshifts of z > 1, a quite 
significant fraction of all star formation occurs in lumi- 
nous, dusty galaxies (as much as 70%; iLe Floc'h et al.l 
[200 5). However, this fraction begins to fall rapidly below 
z = 1 and it is already rela tively low by z = 0.4: Figure 
14 of lLe Floc'h et al.l (|2005[ ) indicates that about 20% of 
star formation at that era occurs in LIRGs, and a con- 
sistent value of 10-20% is given bv lPerez-Gonzalez et al.l 
(|2005f) . It is therefore not unexpected that if a small 
fraction of short GRBs were produced by young (<100 
Myr) stars, we might have observed one in a LIRG — 
though it is likely that several other short GRBs would be 
observed to occur in more "typical" star-forming galax- 
ies first. (Whether this has in fact occurred is a matter 
of debate: see SJTJ) 

The fraction of all stellar mass i n LIR Gs at moderate 
redshift is given by ICaputi et al.l (|2006( ) , who estimate 
that 24% of the stellar mass at z = 0.5-1.0 is in LIRGs 
and ULIRGs. The fraction at lower redshifts is not given. 

While extreme, the qualitative properties of Gl are not com- 
pletely without precedent among short GRB hosts. T he host of 
GRB 070724 (Berger et_ah 2009,; Kocevski et al.|[20lOl ) in partic- 
ular, despite having an SFR and stellar mass an order of mag- 
nitude lower than those of Gl, is also well above the average 
among short-GRB hosts in both these properties. The host of 
GRB 070724 is also significantly dust-extinguished, and in fact has 
a dust-reddened optical afterglow. 



but Figure 8 of ISalim et al.l (|2009f) suggests it is similar 
(within a factor of a few) to that at z w 0.5-1.0; the 
number of LIRGs does not begin to fall precipitously to 
the local density until about z < 0.3, when they essen- 
tially disappear in deep surveys. The fraction of all stars 
cosmologically that are present in LIRGs in this redshift 
range is therefore about 10-20%, quite comparable to the 
cosmological fraction of young stars in LIRGs. Even by 
a posteriori arguments, then, we are unable to associate 
this burst with a specific stellar population age with any 
confidence. 

We can also examine the galaxy's morphological prop- 
erties for clues about the GRB's origin. At z ~ 0, for 
example, LIRGs are predominantly or entirely the prod- 
ucts of major mergers: transient and (relatively) rare 
phenomena in the life of a galaxy with clear observational 
signatures, such as nuclear starbursts and tidal features. 
(This is not the case at z « 1, when mer ging systems 
becom e a minority among LIRGs/ULIRGs: IZheng et al.l 
12004 ) Independent of the SED, then, a GRB from within 
a merger system could be seen as an indication of asso- 
ciation with a progenitor with age similar to the merger 
itself, on the grounds that the explosion of an older star 
by chance during the short-lived merger phase would be 
an unlikely coincidence. No such signatures are present; 
while our resolution is not sufficient to robustly rule out 
a merger remnant, the majority of the galaxy's fiux ap- 
pears to be associated with an ordinary exponential disk. 
The galaxy is slightly asymmetric, in that the region of 
maximum fiux is displaced by ~ 1 kpc from the center of 
the outer isophotes, but this does not appear to reflect 
the detailed properties of the galaxy. Similarly, there are 
no obvious disturbances in the structure of the rotation 
curve (Figure 2]). 

Finally, we can examine the physical characteristics of 
the galaxy from our SED modeling. Evidence for re- 
cent, significant elevation in the SFR above its histori- 
cal average (equivalently, a large specific star-formation 
rate SFR/M* ^ 1/t) would also be indicative that 
the GRB caught this galaxy in an unusual, transient 
part of its history, which is much more likely for a 
young progenitor than an older one. But again, there 
is no such evidence: the current and past SFRs are 
comparable, and the specific SFR is a relatively unre- 
markable (SFR/M* wO.5 Gyr-i). This is the largest 
value for a short GRB host yet (SFR/M ^ ranges be- 
tween 0-0.2 Gyr "^ using the SFRs from iBergeil l2009l 
and masses from iLeibler fc Berged [2010l ). but only by 
a factor of 2-3. The value does overlap the distribution 
of (massive-star associated) long GRB hosts (typic al val- 
ues of SFR/M^, range from 0.1- 50 Gyr— 1, Savaglio ^et al.l 
l2009HCastro Ceron et all|201 01. but only at the low end: 
it is not an extreme value indicating an unambiguously 
star-bursting galaxy. 

All available evidence, therefore, suggests that this is a 
relatively ordinary z « 0.4 galaxy when weighted either 
by star formation or by stellar mass. Furthermore, it 
is (probably) being observed at a fairly typical moment 
in its history, with no clear indication of a large, recent 
increase in the SFR related to a merger or other transient 
event. The properties of the galaxy give no clear means of 
distinguishing a short-lived from a long-lived progenitor 
of this system. 



The LIRG Host of Short GRB 100206A 



13 



5. CONCLUSIONS 

Seen as an isolated event, the age of the progenitor of 
GRB 100206A is frustratingly inconclusive: essentially 
every diagnostic we have applied to the host system is 
consistent with both a short-lived progenitor (<10 Myr) 
and a long-lived one (> 1 Gyr). While the galaxy has 
a large current SFR (~ 30 Mq yr~^), this is matched 
by a substantial pre-existing stellar mass (> 10^^ M0). 
This, and a largely ordinary disk-like morphology, sug- 
gest a steady star-formation history and give no direct 
indication of the progenitor's age. 

Seen in the broader context of other short GRB host- 
galaxy work, however, our conclusions can be seen as gen- 
erally supportive of the prevailing paradigm: as discussed 
in i i4.2[ the significant amount of stellar mass present in 
LIRG systems at z > 0.3 suggests that it should be no 
surprise that a few percent of short GRBs should occur 
within these massive and actively star-forming galaxies, 
even if their progenitor is exclusively a long-lived object 
such as a neutron-star binary. The fraction should be 
even higher if some of the detection biases are consid- 
ered: preferential ejection of the progenitor system from 
low- mass host galaxies (due to supernova "kicks"), as 
well as preferential detection of X-ray or optical after- 
glow from systems in dense environments, would both 
have much less impact on massive galaxies such as Gl 
than on more typical lower-mass galaxies (which would 
be actively selected against), favoring the discovery of 
this type of host. Given what has become a quite large 
short-GRB sample (46 objects with X-ray or optical lo- 
calizations and at least 20 likely host galaxies), the oc- 
currence of an event in a LIRG is not only unsurprising 
for an older progenitor, but could have been predicted. 

This, coupled to the lack of any definitive signature 
tying GRB 100206A to recent star formation in its host 
galaxy despite careful scrutiny of several possible lines 
of argument (including the SED-inferred star-formation 
history and morphology), can be interpreted as evidence 
that this event is the product of a relatively mature pro- 
genitor system in the galaxy. Given the precedent set by 
other short GRBs, the most natural explanation for this 
burst is that it was produced in the same manner as the 
rest of the population — for which a compelling statisti- 
cal case is now building for a generally old (> 100 Myr) 
progenitor, even if the delay-time distribution is indeed 
broad. 

Of course, a younger progenitor is not ruled out. If 
this event were the result of a young progenitor (and 
especially a massive star), it is interesting to note that 
the high metallicity of this system — at least Solar and 
perhaps significantly higher — would nevertheless be in- 
dicative of an origin distinct from that of long GRBs, 
which quite notably avoid dusty, massive LIRG-like sys- 
tems in the local universe. In that case this event would 
be more indicative of a broad delay-time distribution for 
the short-GRB population, rather than an origin from 
the same population that p roduces LGRBs du e to, for 
example, off-axis collapsars (jLazzati et al.|[20To|) . While 
it is possible that a few long-d uration GRBs may oc - 
cur in high-metallicity regions (jLevesgue et al. 2010b'), 
these appear to b e relatively rare f|Modiaz et al. 2008; : 
iStanek et al1[2006h . so even if this event were to be asso- 
ciated with the youngest stellar population in this galaxy. 



it would nevertheless continue to support a distinction 
between the progenitor systems of short and long GRBs. 
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